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Motivated by neutron scattering data, we develop a model of electrons interacting with a magnetic resonance 
and use it to analyze angle resolved photoemission (ARPES) and tunneling data in the superconducting state 
of Bi2Sr2CaCu20g+6. We not only can explain the peak-dip-hump structure observed near the (71", 0) point, 
and its particle-hole asymmetry as seen in SIN tunneling spectra, but also its evolution throughout the 
Brillouin zone, including a velocity 'kink' near the d-wave node. 
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Recent advances in the momentum resolution of 
ARPES spectroscopy have led to a detailed mapping 
of the spectral lineshapc in the high-T c superconductor 
Bi 2 Sr 2 CaCu20 8+ 5 (BSCCO) throughout the Brillouin 
zone. |]|| In contrast to normal state data where well 
defined excitations do not exist, quasiparticle peaks were 
identified below T c along the entire Fermi surface. || 
Moreover, it has been known for some time that near the 
(71-, 0) (M) point of the zone, the spectral function shows 
an anomalous lineshape, the so called 'peak-dip-hump' 
structure. 

§1 The new data indicate a) near the d-wave 
node of the superconducting gap, the dispersion shows a 
characteristic 'kink' feature: for \u>\ < uJkink, the spec- 
tra exhibit sharp peaks with a weaker dispersion; above 
this, broad peaks with a stronger dispersion; b) away 
from the node, the dispersion kink develops into a 'break'; 
the two resulting branches are separated by an energy 
gap, and overlap in momentum space; c) towards M, the 
break evolves into a pronounced spectral 'dip' separating 
the almost dispersionless quasiparticle branch from the 
weakly dispersing high energy branch (the 'hump'); d) 
the kink, break, and dip features all occur at the same 
energy, independent of position in the zone. [[lj 

Features similar to the ARPES spectrum near the M 
point were earlier observed in tunneling spectroscopy. || 
Experimental SIN junctions on BSCCO show a charac- 
teristic asymmetry, with a more pronounced dip-hump 
structure on the occupied side. On the other hand, 
SIS junctions reveal a strong dip-hump feature on both 
bias sides. || 

There have been several theoretical treatments which 
assigned the anomalous ARPES lineshape near the M 
point of the zone to the coupling between spin fluctua- 
tions and electrons. [p|-|l2| Here, we are able to explain 
features a)-d) of the ARPES data, as well as the SIN 
tunneling asymmetry, in terms of the combined effect of 
A) the flat electronic dispersion near the M point of the 
zone and B) coupling of the fcrmionic degrees of freedom 
to a bosonic mode which is sharp in energy and peaked 
in momentum near Q = (tt,tt). Our main result is that 
the anomalous features in the dispersion and lineshape 



for all points in the zone have the same origin. 

A resonance mode with these characteristics is ob- 
served in inelastic neutron scattering experiments in bi- 
layer cuprates in the superconducting state. |l3|,@ The 
neutron resonance lies below a gapped continuum, the 
latter having a signal typically a factor of 30 less than 
the maximum at Q at the mode energy. In order to 
extract the essential physics, we concentrate on the mode 
part and neglect the continuum. The latter contributes 
mainly to additional damping at higher energies. We 
treat the mode in a semi-phenomenological way, taking 
the relevant parameters from experiment. We then calcu- 
late the resulting electronic self energy to second order in 
the coupling constant. From the self energy we directly 
obtain the spectral function measured by ARPES, which 
we then use to calculate the tunneling conductance. 

The retarded self energy on the real energy axis is given 

by H 
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with A* B(k,e) = ZflZo -q,e- u)B(q,u), 

fj Pauli matrices in particle-hole space, and g the effec- 
tive coupling constant. The Keldysh (K) components 
are given in terms of retarded (i?) and advanced (^4) 
functions by G K = (G R - G A )(1 - 2/) and X K = 
(x R ~ X A )(1 + 26), with the usual Bose (b) and Fermi 
(/) distribution functions. 

The model for the mode is based on measurements of 
the spin susceptibility from inelastic neutron scattering 
experiments. fi~I|| The mode energy will be denoted by 
£1 and its energy width is assumed to be irrelevant for 
the self energy. This assumption will be confirmed later. 
This leads to the following model for the mode part of 
the susceptibility 
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Here f(q) describes the momentum dependence of the 
mode and is assumed to be enhanced at the (it, it) point. 
Using the correlation length £ we write it as 
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Experimentally the energy integrated susceptibility at 
the (it, 7r) wavevector, t^Xo-i was determined to be 0.95/j,% 
per plane for BSCCO, Q leading to Xq = 0.3^|>. For 
the correlation length, we take a conservative estimate of 
£ = 2a. This corresponds to a full width half maximum 
of 0.26A -1 , as observed in YBCO, but somewhat smaller 
than that estimated for BSCCO (0.52A" 1 ) 0] which we 
feel is somewhat broad. The mode energy was chosen 
to be SI = 39 meV, which represents the reported values 
between 35 and 43 meV. Q 

Though x is 'renormalized', we use a bare G in Eq. 1. 
This is the same approximation as in Ref. p7y, where it 
was shown that this is better than using renormalized 
G with vertex corrections neglected. This is unlike the 
electron-phonon problem, where Migdal's theorem ap- 
plies. We take the success of explaining the experimental 
features as strong support of this approximation. 

The bare Green's functions with normal state disper- 
sion £e, gap function Ag, and excitation energy E? = 
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where a n - (1 + £%l 0n = (1 - k/ E k)/ 2 > 
Q'12 = —0i2 = ~A-k/2Ej:, etc. For the normal state 
dispersion we use a six-parameter tight binding fit. |Fq| 
We neglect bilayer splitting, as experiments suggest it 
is absent in BSCCO. || A characteristic feature of this 
dispersion is a flat band with a saddlepoint at M with en- 
ergy £m = —34 meV. The superconducting gap function 
is taken to be the d-wave = Ao(cosfca; — cosk y )/2 
with a maximal gap value of A = 35 meV. We have 
chosen r — 5 meV as an intrinsic lifetime broadening. 
The coupling constant relevant for our model is 9 2 Xqi 

chosen to be 0.125 eV 2 /i 2 B . Given a value Xn ~ O-^Ms) 
this corresponds to g — 0.65 eV, the same value as used 
in previous spin fluctuation work. fl9[| We performed the 
w-integration in Eq. 1 analytically and the correlation 
product in momentum space via fast Fourier transform, 
using 256 x 256 points in the Brillouin zone. 

In Fig. 1 we show the renormalization function Z(e) 
1 — Eo(e)/e, where Eo is the to component of the E ma- 
trix. Since the q integral in Eq. 1 is dominated by the 
regions around the M point where the band is flat and 
close to the chemical potential, there are features in the 
imaginary part of the self-energy connected with the two 
extremal energies A and Em 
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tures do not show dispersion, but a change in magnitude 
with position in the zone which is determined by the mo- 
mentum width of the mode. This is the central result of 
this paper. More generally, the imaginary part of the self 
energy is enhanced between the values 
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FIG. 1. Renormalization function Z(e) at the M=(n,0) 
point (solid line) and at the node (0.36-7T, 0.367r) (dashed line) 
for T = 40K. 

t\ = Ao + SI and 62 = Em + ^- For £m approaching 
the chemical potential, Em approaches Ao resulting in a 
peak-like feature in the self energy. For our case, e\ = 74 
meV and £2 = 88 meV. Because the spectral weight of 
the mode is maximal near Q — (tt,tt), the M points of 
the zone, which are connected by Q, benefit mostly. This 
results in stronger features in the self energy near the M 
points compared to e.g. the nodal points. The peaked 
structure in the imaginary part of the self energy results 
(via Kramers-Kronig relations) in an enhancement of the 
real part of the renormalization function for |e| < ex, and 
a reduction of it for |e| > €2, as shown in Fig. 1. This 
leads to a renormalization of the low-energy dispersion 
of the spectra compared to the high energy part. Since 
the experimental energy width of the neutron resonance 
is smaller than the variation in energy of typical features 
in the self energy, this confirms our assumption that the 
energy width of the mode is not relevant. 
The spectral function is obtained by 
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In Fig. 2 we show the spectral functions for momentum 
cuts through the M point towards Y = (tt, it) (MY cut), 
and parallel to MY through the nodal point. Due to 
particle-hole coherence factors, there are quasiparticlc 
peaks at M on both sides of the chemical potential. On 
the negative energy side, the peak is more pronounced 
since £m is negative, and a strong dip feature is present. 
The asymmetry of the dip feature is a combined effect 
of the T3 component of E, which introduces particle- hole 
asymmetry, and the inherent particle-hole asymmetry of 
the band structure near the M point. Going from M 
towards the Fermi surface (Fig. 2, bottom), the hump 
feature quickly loses weig ht as observed in ARPES. @] 
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FIG. 2. Spectral function A(k, e) for T = 40K. Top panel: 
for k points along the cut parallel to MY through the order 
parameter node, from (0.207T, 0.367r) to (0.517T, 0.367r). Bot- 
tom panel: for k points along the MY cut from (0, n) to 
(0.437r,7r). 

In the top panel of Fig. 2 we show spectra for a cut paral- 
lel to MY through the order parameter node at the Fermi 
surface. Near the node there is only one peak crossing 
the chemical potential. The dip-hump features are very 
weak near the node and are presumably overshadowed 
by the additional lifetime effects due to the continuum 
part of the spin susceptibility. Note the much broader 
peaks for higher energies, |e| > 80 meV, compared to the 
sharper peaks near the chemical potential, as observed 
in ARPES experiments. @-[| 

In Fig. 3 we present our results for the dispersion ob- 
tained from the maxima of the occupied part of the spec- 
tral function, A(k,e)f(e). Near the M point we observe 
an almost dispersionless strong peak feature at roughly 
the gap energy — Ao, and a weaker hump feature at 
slightly below —£2, consistent with experimental finding 
c). The peak feature, which without interaction with the 
mode would be at Em, is pushed towards the chemical 
potential, thus ending up close to Ao for not too small 
coupling constants. The position of the hump feature 
is strongly dependent on the coupling constant. We ad- 
justed g to reproduce the experimental value of about 
-130 meV for the hump feature at M; this choice also 
results in the weak dispersion of the hump feature as ob- 
served in experiment. [Q As one goes away from M, the 
dispersion of the hump extents further below —62 and 
the peak starts to show dispersion, until a characteristic 
break in the dispersion with a jump at ^-80 meV de- 
velops. This is exactly the experimental finding in point 
b). Note the stability of the characteristic -80 meV en- 
ergy value for the break/dip feature throughout the zone. 
This is a result of the dominance of the region near M 
in the q sum in Eq. 1, which sets the energy scale. Thus 
we confirm point d) of the experimental findings. 
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FIG. 3. Calculated dispersion at T = 40K obtained from 
the maxima of the occupied part of the spectral function for 
cuts along k x for fixed k y as indicated. Momenta are in units 
of 7T. The size of the symbols represents the peak intensity. 
The coupling constant is g = 0.65 eV except in panel 1, where 
we show for comparison with panel 2 results for g = 0.39 eV. 
Note striking similarity of this plot to the data of Ref . Jlj . 

As the nodal point is approached, the self energy becomes 
weaker due to the momentum dependence of the mode. 
The sudden change in the lincwidth for a cut parallel to 
MY through the node (panel 2 of Fig. 3), as discussed in 
Fig. 2, occurs around -80 meV, in accordance with point 
a). We still observe a weak break feature, which will be 
smeared out by additional lifetime broadening from the 
continuum part of the susceptibility. This weak break 
is also reduced for a smaller coupling (panel 1), or if 
the Lorentzian in Eq. 3 is replaced by a Gaussian. Note 
that in accordance with experiments, the velocity near 
the nodal point is reduced compared to that for higher 
binding energies, causing a velocity 'kink'. 

Knowing the spectral function throughout the zone, 
we are able to calculate the tunneling spectra given a 
tunneling matrix element T^-. From the SIN tunneling 
current /(V") one obtains the differential conductance, 
dl/dV. As usual we neglect the energy dependence of 
the SIN matrix element |A/ £ | 2 = 2eJ2p\Tk P \ 2 A N (p,e), 
where Apj is the spectral function of the normal metal. 
The SIN tunneling current is then given by 

/°° fj f 
A(k,e){f(e)-f(e + eV)} (6) 
-00 zi7r 

k 

In the top panels of Fig. 4, we show results for the SIN 
dl/dV for several coupling strengths. We model the tun- 
neling matrix element for two extreme cases: for incoher- 
ent tunneling we assume a constant |M^| 2 = M5 2 , whereas 
for coherent tunneling we use |M^| 2 = jMi(cosk x — 
cosk v ) 2 . H 
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FIG. 4. Differential tunneling conductance for SIN (top) 
and SIS (bottom) tunnel junctions for T = 40K. Units are 
eM, 2 for SIN and 2e 2 T? for SIS. Left for isotropic tunnel- 
ing (i=0), right for coherent tunneling (i=l). Curves are for 
3=0.39 eV (dotted), 0.65 eV (full line), and 0.90 eV (dashed). 



In conclusion, we have shown that the momentum dis- 
persion of the ARPES spectra, as detailed in recent ex- 
periments, can be explained by a simple model which has 
as components A) a flat band region near the chemical 
potential in the normal state dispersion near the (n, 0) 
point of the zone; B) a nearly dispersionless bosonic mode 
which is peaked in momentum near the (jr, n) point, and 
which interacts with the fermionic degrees of freedom. 
The theoretical tunneling spectra obtained with the same 
parameter set are consistent with the experimental find- 
ings of an asymmetry of the peak-dip-hump structure in 
SIN tunneling spectra. 
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In both cases, we observe a clear asymmetry with a dip- 
hump structure on the negative bias side and a very weak 
feature on the positive side of the spectrum, as in exper- 
iment 0. The pronounced asymmetry is a result of the 
shallow band near the M point, £m ~ — which en- 
hances the coupling to the resonance mode for populated 
states. Note that the hump position is strongly depen- 
dent on the coupling constant in contrast to the position 
of the dip minimum. The asymmetry in the peak height 
on either side of the spectrum is sensitive to the coupling 
constant too. In weak coupling the negative bias peak is 
higher due to the Van Hove singularity at the M point. 
For stronger coupling the pronounced dip at negative bias 
reduces the height of the coherence peak on this side and 
shifts the hump to higher energies. For g=0.65 eV (full 
lines in Fig. 4) the peaks at positive and negative bias 
have roughly the same height, as in experiment [jTj. 

For an SIS junction, the single particle tunneling cur- 
rent is given in terms of the spectral functions by 

/°° df 
A(k,e)A(p,e + eV) 
-OO ^ 

kp 

x{f(e)-f(e + eV)} (7) 

Again we show results for incoherent tunneling (|Tg^| = 
Tq) and for coherent tunneling with conserved parallel 
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momentum, |T^| 2 = ^jTf(cos k x — cos k y ) i 5(k^ — 
Q We show the SIS dl/dV in the bottom panels of 
Fig. 4. Our theoretical SIS curves for incoherent tun- 
neling resemble very closely the experimental results for 
BSCCO, H unlike for coherent tunneling which exhibits 
negative dl/dV regions due to the strong anisotropy of 
Tjt-. Note that the dip-hump feature is strong on both 
sides for an SIS junction in contrast to the SIN results. 
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